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Abstract 



Cumulant oscillations, or Hq inoinent oscillations, appear if the KNO multipar- 

Q^' tide distribution decreases at large z, z = n/ < n >^ faster than the exponential, 

' , exp{—Dz^), with fi > 1. In nucleus-nucleus interactions this behaviour is related 

Qj I to the limitation in the average number of elementary central collisions (or aver- 

1-^ ' age number of strings centrally produced), due to the finite number of nucleons 

involved. Colour deconfinement, via percolating string fusion, will drastically de- 

K^ \ crease the fraction of centrally produced strings and increase the cut-off parameter 

;_( ■ //: Moment oscillations will be displaced to smaller q and the width of the KNO 



distribution and forward-backward particle correlations will become smaller. 
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In the framework of perturbative QCD it was theoretically predicted, sometime ago, 
that the factorial cumulants Kq (or, equivalently, the moments Hq = Kq/ Fq, Fq being 
the factorial moments) of the multiparticle distribution should present oscillations in 
sign as a function of g [|l|. The predictions turned out to be confirmed by experiment in 
e+e~ annihilations 0, in pp interactions |]^ and in hadron-nucleus and nucleus- nucleus 
collisions 0. It is now clear that such oscillations are present in all known high energy 
processes. 

Recently, in |p, it was shown that a necessary condition for a particle distribution 
P(n, < 77, >), < n > being the average multiplicity, with positive two particle correlation, 
K2 > 0, to have oscillations in Kq is to exist, asymptotically, a KNO distribution "^{z), 
M , of the form: 



< n > P{n, < n>) > ^(2;) > exp[-Dz'^] , (1) 

<n>^oo ^>1 

z=n/ <n>=const. 

z = n/ < n > being the scaling variable, D and /i positive parameters, with 

fi>l. (2) 

Most of the existing popular parametrizations (of the Negative Binomial Distibution 
family) do not satisfy (^ and (||). They behave exponentially at large z, fi = 1, and, 
not surprisingly, they do not show oscillations in q: Kq > for all values of q. 

The importance of the large z behaviour of the distributions to generate the oscilla- 
tions can be easily seen from the expression relating Kq to Fq {Fq =< n{n — l)....{n — 
g + 1) > / < n >'?)) and Kp, p < q: 

K, = Fq- J2 nU...^n^..-, (3) 

ri,r2,... 

with ri + r2 + ... = q, ri > 1, rj+i > r^ and ri < q. The -P^i r2 ... ^^^ ^^^ positive 
combinatorial factors associated to the partition of an integer q into integers ri, r2, ... . 
The meaning of (|^) is straightforward. To obtain the q particle cumulant Kq one has to 
substract from the q particle factorial moment Fq (integrated q particle inclusive density) 
all the q — 1, q — 2,... particle cumulants (integrated inclusive particle correlations) in 
all clustering combinations. We have, in particular (note that we are using normalized 
moments, Ki = Fi = 1) , 

K2 = F2-1 

Ks = F,- 3K2 - 1 (4) 

Ki = Fi- AKs - 3Ki -QK2-I 

7^5 = F5 - 57^4 - IOK3K2 - IOK3 - 15K^ - IOK2 - 1. 



If the distribution if of exponential type the Fq grow very fast with g, essentially as 
g!, and all the terms in (^) and (^ remain positive {K2 is assumed positive). If at large 
z the distribution decreases faster than the exponential, /i > 1, the Fq grow slower and 
at some value of g, Kq start becoming negative. But this negative Kq enter in the future 
steps, g + 1, g + 2,. . . giving positive contibutions to the right hand side of (^ and (|). 
At some stage the Kq become again positive, etc.: the oscillations start. 

It can be argued that the required supression at large z is not really physical, in the 
sense that at finite energy and limited experimental acceptance, distributions are always 
cut at large z 0. In e^e~ annihilations there are theoretical reasons to belive that the 
oscillations are physical |1|. We shall argue here that in nucleus- nucleus collisions, at 
least, the oscillations have also a physical origin. 

Succesful models that attempt to explain hadron-hadron, h-h, hadron-nucleus, h- 
A, and nucleus-nucleus, A-B, collisions, making use, or not, of basic information from 
e^e~ annihilations, are multiple scattering models. We shall take as reference the Dual 
Parton Model (DPM), p]], but most of our results do not depend on detailed features of 
a particular model. 

In any multiple scattering model the distribution P{n) of produced particles is the 
result of the superposition of the contributions from elementary inelastic collisions. At 
each elementary collision particles (via string formation, for instance) are emitted with 
a given distribution. There is a certain probability of v elementary collisions to occur. 
In general, we can write 

Pi^) = ^ X^ V{v)p{ni)p{n2)...p{n^) (5) 

with rti + n2 + ■ ■ ■ + Ui, = n, p{ni) being the particle distribution from the i-th el- 
ementary collision, v^(z^) the probability distribution for z/ elementary collisions. The 
parameter u represents, as well, the number of intermediate produced objects: pairs of 
strings, in DPM. In (^ we have assumed that the formed strings emit independently 
and that fluctuactions in the size of the strings are negligeable, or can be reabsorbed in 
the distribution ip{i'). 

Let us now introduce the generating function G{z), 



G{z)=J2i'^ + zrP{n) (6) 

n=0 

with G(0) = 1, such that 

1 cf^Giz 



F = '. 

'*' < n >9 dz'i 

and 

1 dnogG{z) 



^ <n>i dzi 



(7) 

z=0 



(8) 

z=0 



By combinig (|) and (|^) we obtain 

Giz) = J2^i'^)9izr, (9) 



u=l 



where g{z) is the generating functioin for the elementary process. Knowing the elemen- 
tary generating function g{z) (from e~^e~ , string model, Poisson approximation, etc) and 
the elementary collision distribution ip{i') (from multiple scattering combinatories, im- 
pact parameter integrations, etc) the full generating function G{z) can be constructed 
and the moments computed, (^, (0) and (|^). A summary of results is contained in Table 
1. 

For the average multiplicity < n > and the normalized KNO dispersion D/ < n >, 
where D^ =< n^ > — < n >^ we have, 

< n >=< u > n (10) 

and 



2 2 + ^ ' (11) 

where n and d are the average multiplicity and the dispersion of the elementary process, 
respectively. 

Concerning relation (|ll|) an important observation can be made. If fuctuactions in 
the effective number of strings were negligeable, i.e., {< v"^ > — < v >^)/ < v >^— 0, 
the second term in the right hand side of (0) should dominate. The KNO dispersion, 
as < z/ >> 1, should then be smaller than the normalized dispersion of the elementary 
process (say, e~^e~ or Poisson distribution). This is against experiment: D"^/ < n >'^ 
increases with the complexity of the systems involved, from 0.09 in e~^e~ annihilations 
to 0.25 - 0.30 in pp collisions |T^, to larger values in h-A processes, and to 0.8 - 1 



in A-B collisions UTTj |l^ . The conclusion is that in ( [Til) the first term in the right hand 
side, counting the fluctuations in the effective number of strings, is dominant. This is 
particularly true in nucleus-nucleus collisions where, for large A and B, < u >~ 10^ — 10^, 
thus making the second term in the right hand side of ( |TI| ) completely negligeable. What 
we have shown for D/ < n > applies to all the cumulants Kg (see Table 1). At least 
in nucleus- nucleus collisions the KNO particle ditribution function < n > P{n,< n >) 
must be very close to the KNO string distribution function < u > V5(z/, < u >), with 
scaling variable z = n/ < n >':^ v/ < v >. 

It should by now be clear that oscillations in the factorial cumulants Kq (or in the Hq 
moments), at least in nucleus- nucleus collisions, have nothing to do with the elementary 
interactions and cannot be related in a simple way to some perturbative QCD calculations 
(which may apply to the elementary process). We shall try next to explain the origin of 
the oscillations in A-B collisions. 



In nucleus-nucleus A-B collisions, if instead of measuring the unconstrained multiplic- 
ity distributions P{n) one measures the distribution at a fixed impact parameter b (in 
particular in very central collisions, at 6 = 0) the obtained distribution is totaly differ- 
ent. While in the unconstrained situation the KNO distribution is wide {D/ < n >~ 1) 
and roughly independent of the nuclei, the central KNO distribution is very narrow 
{D/ < n ><^ 1) and strongly dependent on the nuclei envolved, specially on the lightest 
nucleus |TT|. These differences can be easily understood from (|TT]). 



If the impact parameter is fixed the fiuctuation in the number of strings is drastically 
reduced, {< u"^ > — < u >^)/ < i^ >^^ 0, and the second term in the right hand side of 



(11) dominates: 



^2 



< n>' 



central 



< u >rn 



2 ' 



(12) 



where < z/ >c is the average number of central elementary collisions. One sees that the 
KNO distribution has to be very narrow {d/n is almost Poisson-like and < v >c may 
be very large), and the width depends on < z/ >c. The value of < z/ >c increases with 
the atomic weight of the nuclei, being limited by the atomic weight of lightest nucleus. 
Independently of A and B one roughly has [ITT] 



< V > / <u >c~ 1/4 



(13) 



A good way of parameterising multiplicity distributions, from e+e^ annihilations to 
A-B collisions, is by using the generalized gamma function for the asymptotic KNO 
function, 



^{z) = 



f^ 



T(k) 



r(/t + i//i) 



K/i 



;"'' ^exp 






(14) 



with 



< z'^ >= ^^"^'i^h i^r\ 



(15) 



and, in particular, < z^ >=< z^ >= 1. The parameters k and fi have to be fixed. 

We shall now impose the constraint that the large n behaviour of the central collision 
particle distribution must be the same as the behaviour of the unconstrained distribution: 



P(n, < n >' 



Pin, < n> 



central 



(16) 



If the central distribution is also parametrized as (p^, with parameters fie and k^, 



it is clear, from (IT3), that fie = yU- However, Kc is very different from k, as the central 



distribution is very narrow. This means, see (|T5|) for q = 2, < z"^ >~ 1 or /tc ^ 1- We 



(17) 



can finally impose (|l^) in a stronger way, see (|T^ 

r(K + l//i) _ <jy> 

In order to have estimates for k and fi we introduce the additional relations, see (|15D, 

T{K + 2/fi) T{k) D' 

T{K + l/fi)T{K+l/li) <r2>2^ ' ^^ 

and make use of experimental information on D/ < n >. Reasonable values for k and /i 
in agreement with [|l^ and [|^, < u > / < u >c— 1/4 and D/ < n >~ 0.9, are: 

K ~ 0.1 
and 

/i ~ 5 (19) 

A few remarks can now be made: 

1) As, from (|19]), /t/x < 1, the KNO function does not turn to zero as 2; — > 0, in 
agreement with experimental data on multiplicity and transverse energy distributions. 



mm 

2) The parameter /i is very large, which means that the particle distribution is 
drasticaly cut at large n, in agreement with multiple scattering models ||ll], |l^, |13|, |14 
and data; 

3) As ji is large, oscillations occur in the cumulants Kq. 
In Fig.l we show a plot of the KNO distribution which agrees with experimental data. 
In Fig. 2 we present Hq as a function of q. In both cases k and /i were fixed at the values 



Let us next suppose that the road to the formation of the quark-gluon plasma in 
nucleus-nucleus collisions is by string fusion ||rB[|, mainly occuring in central collisions, 
where the density of strings is higher. The net result is that the large n tail of the 
multiplicity distribution is further cut at large n, or, in other words, fi is further increased. 
This can directly be seen in ([T7|): an increase of < z/ > / < i/ >c, due to string fusion, is 
translated in an increase in /i. 

In Fig.l and 2 we also show the KNO multiplicity distribution and moments Hq as 
function of q when the ratio < z^ > / < z/ >c is increased from roughly 1/4 to 1/2. The 
oscillations tend to start earlier. 
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The same kind of effect occurs if one is triggering on a heavy particle hke the J/^. 
The n-tail of the multiphcity distribution is further cut and therefore the oscillations tend 
to start earlier, as it is seen in Fig. 2 where it is ploted Hq for the multiplicity distribution 
in S-U collisions at y/s = 19.4 GeV when a J/\E' is triggered. Also, It is shown the result 
when Drell-Yan pairs are triggered instead of J/\E'. 

We look now to the forward-backward correlations. If one fixes the number of particles 
forward, and studies the distribution in the backward hemisphere, one approximately has 



with 



< ub >= a + bnp 



dIb/dI 



(20) 
(21) 






--< ubUf > — < Ub >< rip > , (22) 

and Dp being the variance in the forward hemisphere. Within the spirit of our approx- 
imation of keeping only fluctuations in the number of elementary collisions, we assume 
that particles are produced symmetrically in center of mass rapidity, see |TH] , to obtain 



and, see (11), 






< Z/^ > - < Z/ >' 

< z/ >2 



< Ub >< i^F > 



(23) 



< z/^ > 



<u >' 



(p 



D'p = + — < n^ >< n^ > , 24 

where the factor 2 in the second term in the right hand side of (p^) corresponds to 
particles being independently emitted in rapidity in each elementary collision. 

It is clear that in an unconstrained nucleus- nucleus, A-A, collision, from (pOD, ( pT]) 
and (p2D , the forward-backward correlation parameter b is essentially 1. If we trigger in 
a central collision, <z/^> — <z/>^/<z/ >^^ 0, as seen before, and the parameter b 



drops to zero. For instance, in the case of percolating fusion |jT^, as < z/ > / < z/ >t 
1/2, if one triggers on forward events with 



nF>2<nF> 



(25) 
but as we have 



one should already obtain 6 ~ 0. This agrees with the results of [| 
here strong fusion the net effect is more spectacular. 
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Figure captions 



Fig.l. KNO negative particle distributions in heavy nucleus-nucleus collisions. Contin- 
uous curve: present energy situation, k = 0.1, /i = 5 in eq. (p!4|). Dashed curve: the 



same distribution for percolating string fusion, with k = 0.1, /i = 10, corresponding to 
about 50% reduction in the density of central strings. 



Fig. 2. The moments Hg = Kg/Fg for g = 2, 3, 4, 5, 6, 7 in the case of nucleus-nucleus col- 
lisions without fusion (/i = 5, black circles), and with fusion {fi = 10, open circles). Also 
represented with crosses (squares) the Hg deduced from the multiplicity distributions of 
S-U collisions at ^/s = 19.4 GeV when a J/\l/ (a Drell-Yan pair) is triggered. 
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Table 1: The factorial cumulants Kq for g = 2, 3 and 4, obtained from 



and 



When the fluctuations in the number of strings are negligeable the cumulants are the 
cumulants resulting from < u > independent sources: Kq = Kq/ < v >'^~^. When string 
fluctuations dominate, as in heavy nuclei collisions, only the terms inside square brackets 
are important. In that limit, the factorial cumulants Kq are given by the cumulants of 
the string distribution. In that same limit the factorial moments Fq are given by the 
moments < v'^ > / < v >^ oi the string distribution. For definition of moments see, for 
instance, 1116 
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